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Abstract 
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I. INTRODUCTION 

Sampling complex integral weights of the form e~ s is important for many applications 
in high energy physics. Complex Langevin equations have been proposed for this purpose 
[1, 2]. Many theoretical questions are still to be answered for complex Langevin equations. 
Most importantly, it is not known a priori if a complex Langevin simulation will converge 
to a stationary distribution. Even if it does, there are multiple stationary distributions that 
satisfy stationary limit equations (Fokker- Planck equation [3] or Schwinger-Dyson equation 
[4]) and it is not known how to identify to which one of these that the complex Langevin 
equation converges. See [4] for a discussion of these issues and references to literature. 

Identification of the particular stationary distribution sampled by the complex Langevin 
equation needs input from the dynamics of the complex Langevin simulation. Analysis 
depending only on the quantities in the stationary limit has not been successful in differen- 
tiating between stationary distributions. To account for this fact, we introduce an effective 
potential approach to complex Langevin equations. Even though this effective potential 
does not have all the properties of the conventional quantum field theory effective action, 
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we show that it governs the probability that the spacetime and Langevin time average of 
the fields takes on specific values. In the infinite Langevin time limit, it is shown that the 
spacetime and Langevin time average of the field must be equal to the minimum of the 
effective potential. As in quantum field theory, a loop expansion for this effective action 
can be defined. The parameter that defines loops is a factor that multiplies the noise in 
the complex Langevin equation. The same parameter is shown to count derivatives in the 
Schwinger-Dyson equation derived from the action S. This point is important, because it is 
known that the stationary distributions of the complex Langevin equation are complexified 
path integral solutions of the Schwinger-Dyson equation, i.e. see [4]. The main idea of this 
paper is to use the minimum of the effective action to identify the stationary distribution 
to which the complex Langevin equation converges by matching the loop expansion with 
the derivative expansion. This provides an approach to addressing the multiple stationary 
distributions problem. 

In two recent studies [5, 6], the classical fixed point structure of complex Langevin equa- 
tions were shown to be important for the convergence behavior. Our results will also parallel 
this result since to leading order the minima of the effective potential will be given by classi- 
cal fixed points. The next order will be used to remove the degeneracy and this will lead us 
to the identification of the stationary distribution that is sampled by the complex Langevin 
simulation. 

We start by reviewing relevant results in Section II. We then introduce the effective 
potential in a lattice setting, discuss its properties and calculate it to one loop order (Section 
III). Before we conclude, we give a simple application of these ideas (Section IV). 

One point we omit in this paper is the discussion of the continuum limit. The question 
of renormalizability is a nontrivial problem and deserves a separate study of its own. 

II. STATIONARY DISTRIBUTIONS OF COMPLEX LANGEVIN EQUATIONS 

In this section we review the relevant results of [4]. Even though, our ultimate goal in 
this paper is to do calculations in a lattice setting, the arguments in this section are given 
in the continuum for ease of notation. The lattice versions of these results are either in [4] 
or [7]. 
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Given a complex action S of a scalar field (f>(x) in a Euclidean space, 
S[<f>] = J d d x l -4>{x) (-□ + m 2 ) <t>{x) + jfa) 



we want to study the stochastic dynamics given by 

(x, r) 8S[4>] 



(1) 



where r/ R and rj 1 are independent Gaussian noises with correlators given by 

( v r (x,t)) = 0, 

(r] R (x, r) r] R (x', r')> = 2n a 5 d (x - x') 5(r - r'), (3) 

and 

(r ] I (x,r))=0, 

{rf{x, r) rf{?t, r')> =2Q(a- 1) 5 d (x - x') 5{r - r'). (4) 

a and are real constants and a > 1. = 1 corresponds to the complex Langevin 
equation. We will not set f2 to any particular value for the moment. Q will define orders in 
the approximations we do below. Equation (2) can be written as a stochastic system of two 
real fields, <f>(x, r) = <p R {x, r) + i0 7 (x, r). In fact, this is what we will do in the rest of this 
paper. 

Note that most of the time in the literature a is set to 1. There is only a real noise in the 
complex Langevin equation. As shown in [4], this may lead to initial condition dependence 
of the stationary distribution sampled by the complex Langevin equation. This happens 
because for some actions the complex Langevin dynamics may be confined to sampling 
in lower or upper half planes depending on the initial condition, not having access to all 
of the complex plane. In this paper, we want to avoid such initial condition dependence 
and assume that the stationary distribution sampled by the complex Langevin equation is 
unique. Furthermore, setting a to some value larger then 1 will be important for deriving a 
path integral representation of the generating functional for the complex Langevin process 
in the next section. 

Given some initial conditions for the field, Langevin equation (2) has a time dependent 
probability distribution, see e.g [8], 

P{t^{x)^{x)) = (8(^{x,t)-^ (5) 
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which can be shown to satisfy a Fokker-Planck equation, 



dP (r, 0« 0') 



dr 



+ 



( ? 5P Rc 5S ^ P 

8(j) R (x) 5(j)(x) 

n { a-i)J^ + iJMp 



(6) 



8(f) 1 {x) \ v 7 5(f) 1 (x) ' ""84>{x) 
The existence and uniqueness of stationary distributions for this Fokker-Planck equation is 
still an unanswered question. 

We continue by defining a generating functional for equal time correlation functions 



(7) 

where J[x) is a real source field. Here the brackets denote an average over the noise. Field 
variables are considered as functionals of the noise through equation (2), hence we include a 
superscript r\. The effect of noise appears in the stochastic average through the probability 
distribution P (see e.g. [8]), 

Z[J, T )= J [d<p R ] [#'] P(T, (f) R , (f) 1 ) e Sd d ^{x,r)+i4>'{x,r))j{x)_ (g) 

This is now an integral over all field configurations weighted by P. This point will be treated 
in more detail in the next section. Using the Fokker-Planck equation (6), one can show that 
the following equation holds 



dZ[J,r) 
dr 



flJ{x) - 



&S[<j>] ( 5 



Z[J,t). 



(9) 



6<f>(x) \5J 

If there is a stationary distribution, the left hand side is zero. When = 1, functional 
derivatives of the right hand side with respect to J produces Schwinger-Dyson identities for 
the action S[(f>], since the integrand must vanish due to the arbitrariness of J. Q plays the 
role of h in quantum field theory. It is easy to see this by defining a making a change of 
variables J(x) — > J(x)/Q} Then, the stationary distribution condition becomes 



d d x 



j(x) - 5 -^& (n— 



Z[J,r) = 0. 



(10) 



Q counts the number of functional derivatives, as h would. As is well known, see [9], a 
derivative expansion of Schwinger-Dyson equations is equivalent to a loop expansion of the 
associated path integral around a local minimum. 



1 In terms of this new variable, the correlation functions are given by functional derivatives of Z[J,t) with 
respect to Qjj- 
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If one considers zero dimensional theories, an explicit construction of the stationary equal 
time correlation functional can be given. Generalization to a lattice is straightforward, see 
[4]. Now, the action S is a function of 0. Equation (9) becomes 

^-nfZ U , T) -M(l)zU,r). (11) 



dr ' d(j) \dj 

In particular, the generating function for the stationary distribution satisfies 

dS f d 



b^\4) ZM) = q1ZJ}) - (12) 

We can solve this equation following [10]. Our ansatz is 

Zr(j)= Jd4>G(4>)j+, (13) 
where T is a contour over the complex plane. Inserting this into equation (12) one gets: 



= -QG(<f>)e>+\ gr + J- 
This equation can be solved for 



e j(p . (14) 



G(0)=e"^, (15) 

and r is contour that connects the zeros of e~ s ^ + ^ on the complex plane. As advertised, £1 
plays the role of h. Since equation (12) is linear, any linear combination of the (independent) 
solutions will also be a solution, 

ZrfO') = £ ar ' Zr 'k')' (16) 

where Tj define an independent subset of paths T. 2 The only constraint on the coefficients 
at this point is the normalization condition for correct distribution interpretation, namely 
that Z st (0) = 1. We restricted the form of the generating function, but did not identify it 
completely. Examples are given in [4]. 

In this paper, we use an effective action approach to identify the stationary distribution. 
The Q expansion gives a small noise expansion for the complex Langevin equation and a 
loop expansion of the stationary distribution generating function. Matching of both series 
order by order will tell us the exact form of the stationary distribution. 



Essentially the same result is derived by another method in [3]. 
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III. EFFECTIVE ACTION FOR COMPLEX LANGEVIN EQUATION 



In this section, we will define the effective action for complex Langevin equation and do 
an Q expansion to gather information about the stationary distribution. We choose to use 
a lattice regularization. We work with an Ito type discretization for equation (2). 

First we define our notation. We introduce a Euclidean spacetime lattice with isotropic 
lattice spacing a. Lattice points are denoted by a d dimensional vector n = (ni, . . . , n^), 
where n» = 0, . . . , (N — 1). The coordinate of a lattice point is given by Our 
convention is to denote Euclidean time with xj. Periodic boundary conditions are assumed 
for the space-time lattice. We discretize Langevin time r with At intervals. Langevin 
simulation is assumed to start at r = Tj and run for N T time steps untill Tf = + N t At. 
In terms of dimensionless lattice variables 0„(/t) = a( d ~ 2 )/ 2 0(na, kAt), e = Ar/a 2 , m = am, 
gi = a (J ( 1 -V 2 )+' g u fjR, 1 — y/aAArr]^ 1 , where the noise correlators are now 

(#[(«) *&(«')) = 2Q (a - 1) <W<W, (17) 
the discretized equation (2) in Ito calculus reads 

4> n {K + 1) = 4> n ( K ) - 6 (-□ + m 2 ) 4> n {K) - 6 Mn 1 ^) + ^ + *^n(«)) ■ ( 18 ) 

I 

In the rest of this paper we assume a > 1. To avoid confusions, we note that 0„(/t) is the 
lattice variable at the spacetime point na and Langevin time Tj+/tAr.Here, the dimensionless 
lattice Laplacian is defined by 

d 

(«), (19) 

m k=l 

where e M is a unit vector pointing along the /c-direction. We assume an initial condition is 
set at r = Ti by specifying the field on the spacetime lattice. 

A. Generating Functional for Dynamic Correlation Functions 

First step in our anaylsis is to define a generating functional for dynamic correlation 
functions, 

Zd[J R , J 1 ] = ^exp 
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(20) 



4> v is a functional of the noise fields through the Langevin equation (18). The sum in k 
parameter starts from 1, since we assume an initial condition </> n (0) is given. Unless explicitly 
stated, k sums and products run from to N T — 1. We use subscript d to differentiate Z d 
from the generating functional of equal time correlation functions defined by equation (7). 
Zd is more general than Z since it generates correlation functions of the field variable at 
different times as well as equal times. We will define an effective action using this generating 
functional. Before doing that, we derive a path integral representation for Z&. Although we 
will not follow their work directly, see [11] for a discussion of path integral representations 
of generating functionals based on Ito discretization and [12] for continuum formulation. By 
definition, Z^ is an average over the probability density functions of the Gaussian random 
variables, 



C is a proper normalization constant for each equation it appears in. The next step is 
to rewrite this expression as an integral over field configurations. Using equation (18), 
we can change integration variables from the Gaussian random variables ^{k) and fj^(K) 



Z d [J R , J 1 ] = C J n^(«)d^(«) exp - Y, (^£7 + 

n,K L n , K 




) 



(21) 




Note that due to our discretization, the Jacobian of this 



transformation is a constant. Defining 
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F*(k) = F R (#?(«),#(«)) = Re m 2 UK) +X>#T 1 (k) 



i 



Fi(n) = F 1 (#?(«),&(«)) = Im m 2 n ( K ) +X>#T 1 (*) 



(22) 
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the generating functional Zd becomes 



Z d [j\ J 1 } =C J [#] exp E { J n(^n («) + ^ («)#(«)) 

' _ ?1 K=l 



x exp 



where the measure is 



-Y 



[©#(«) + eF*{K) 



4ae 



+ 



4(a- l)e 



(23) 



[#]=nn*w*w (24) 

There is no integration over </> n (0), since we assume it to be set by initial conditions. The 
exponent in the second line, with the normalization factor C, is the probability density of 
a particular realization of the Langevin process given some initial <f> n (0). We will denote it 
by T (0(iV T )|0(O)). The previously defined probability distribution (5) and T (j>(N T )\<f>(of) 
are very closely related, see e.g. [13]. For T, 1/Q is an overall factor for the exponent, which 
will allow us to define a loop expansion for this generating functional. 

For later convenience, we define the quantity in the exponent of the second line of equation 
(23) as 



4ae 



+ 



1 2 



4(a - l)e 



(25) 



Forgetting the complex Langevin equation roots for one moment, one can think of equation 
(23) as a the generating functional of a field theory on a d + 1 dimensional lattice with an 
action given by equation (25). This is not a usual field theory, because it does not have 
translational invariance due to k = terms appearing in equation (25). One end of the 
lattice (k = 0) is fixed, whereas the other end (k = N T ) is free. This picture will be useful 
later when we discuss the effective potential. 

Before continuing, we note that having a > 1 was a great convenience for this calculation. 
Otherwise, we would start with averages over only the real noise variables and this would 
lead to a difficulty when changing variables in the path integral from noise to field variables. 
The number of field variables 0^' 7 (k) would be half of t) r (k). Therefore, one would need to 
choose half of the field variables as independent ones and express the remaining in terms of 
the independent ones. Having a > 1 saved us from these difficulties. 



9 



B. Effective Action, Effective Potential and Their Interpretation 



Now we do the standart quantum field theoretic procedure to define an effective action 
for this generating functional, 



W[J* J 1 } = In Z d [JK, J 1 }, 

r[0 fi ,0 7 ] = -wijK'iJ 1 '*] + 

where 



N r 



n k=1 



(26) 



#* j (k)=#(k) if JJ(«) = 

0n' J («) = <fe if ■#(«) = #*(*)■ (27) 



By definition, it follows that 



(28) 



For field configurations that are homogenous and static in Langevin time, the effective action 
reduces to "effective potential" , 

vM = &. (29) 

In the next section, we will calculate the effective potential to one loop order (first order in 
Q). 

We note that the effective action defined above does not have all the properties one re- 
members from quantum field theory. In particular, it does not have the energy interpretation 
[14] because the complex Langevin equation (2) as a system of two real Langevin equations 
is not derivable from an action. 3 See [15] for a detailed discussion. However, one useful 
property is still valid. The effective potential governs the probability that the spacetime 
and Langevin time average of the field takes specific values. To show this, we follow closely 
the argument of [15]. 



3 A Langevin system that is derived from an action will have component equations q-i = —M: + f]i- 
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We ask the probability of spacetime and Langevin time average of the field variables 
taking a specific value, 



Prob 



/ 



^ ££**(«) 

n k=1 

[4] t (0^) i0(o)) 5 ( E - ) 5 ( E E - 




T(0(iV T ) |0(O)) JdXWexp 



n k=1 / \ n k=1 



x exp 



27rAR EE^)-^^ 

\ n k=1 

^MEE^^)-^^ 

= J dXW Z d [J R = i2n\ R , J R = z2ttA 7 ] exp [-i2irN d N T (\ R R + A 7 7 )] . (30) 

At this point we will make an approximation, though we will argue that our approximation 
becomes exact in the infinite Langevin time limit. We are looking for a simplification of the 
effective action evaluated in the presence of a constant source field. For this purpose, we 
consider a derivative expansion of the effective action. The first term will be the effective 
potential, as defined above. Other terms will have factors of (lattice) derivatives of the field 
variable (j) n {n). Note that we are now regarding the field variable as the expactation value of 
a "quantum" operator of the lattice field theory given by (23). Consider complex Langevin 
initial conditions that have no space-time dependence. Then, for a constant source field, one 
expects the derivatives along the space-time directions to vanish due to periodic boundary 
conditions and translation invariance. On the Langevin time direction translation invariance 
breaks, because of the existence of an initial condition. Even though this is the case, we will 
argue that in the long Langevin time limit, one can ignore higher order terms and safely 
approximate the effective action by the effective potential. Consider for a moment a complex 
Langevin simulation starting with space-time independent initial conditions. Assuming con- 
vergence to a stationary distribution, in the long Langevin time limit the field expectation 
values will "forget" the initials conditions and approach a stationary value. This means that 
the Langevin time derivatives of the field expectation values will be small, in fact will vanish 
in the long time limit. By definition, one would be able to calculate this effect using the 
generating functional given by (23) by setting the source fields to zero. We are interested in 
constant source fields that are not vanishing. Now let's go back to the Euclidean lattice field 
theory picture defined by (23). For zero source fields, in the very large k region the lattice 
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looks pretty much translationally invariant as argued above. Adding a constant source field 
will not break this invariance. Therefore, vanishing of the Langevin time derivatives must 
still hold. We argue that if this is the case, the effective potential will dominate over all 
higher order terms in the derivative expansion. To motivate this last comment a little more 
rigorously, consider a constant 7 (like the effective potential) and a function (3{k) (like the 
higher order terms in the derivative expansion) that dies as k — > 00. Now consider the 
quantity 

££i0(") 



lim 

7V T — >oo 



(31) 



We will show that this limit is zero, which will prove that our comment is valid, that the 
effective potential dominates over higher terms in the derivative expansion. To prove that 
it is zero, we will show that it is smaller than any positive £. We will need to recall the 
definition of a convergent sequence. A sequence Q K converges to the limit Q if, for any e > 0, 
there exists a P such that \Q K — Q\ < e for k > P. Now we apply this definition to sequence 
P(k), such that e = I7I £. Now there exists a P^ such that \P(k)\ < \^\ £ for k > P$. Then 



lim 

N T ^QO 





lim 1 




N T ^OC I 



< lim — 

N T ~*oo N T 



E 



(5{K) 



+ 



E 

K=P i + l 



(5{k) 



E«=i £(«) 



lim 



7 




lim jl v m 

K = P i +l ' 



< 



lim 



K = P e + l 



7£ 



7 



for any £ > 0. 



(32) 



Therefore, we can safely assume that a constant source field implies constant 4>n( K )- As 
discussed above, this constant 4> J n is the expectation value of the field variable in the long 
Langevin time limit in the presence of a constant source field. Then 

Z d [J R = i2n\ R , J R = i27i\ I ] = exp jiV d iV T ^2vrA R R ' A + i27r\ I ^ I ' x - V[0 R '\ /,A ]) } . 

(33) 

Here the superscpript A is used to denote that fields are functions of the constant source 
field. We will look at infinite Langevin time limit. Then using the method of stationary 
phase 



Prob 



JJZN E E = ^ ) « ex P [-N d N T V(r, + O (1)] . (34) 
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Thus, the effective potential governs the probability distribution of the spacetime and 
Langevin time average of the field. As the simulation time N T goes to infinity, the space- 
time and Langevin time average of the fields must equal to the minimum of the effective 
potential 4 . Our argument for effective potential dominance assumed space-time indenpen- 
dent initial conditions. However, the stationary distribution that the complex Langevin 
equation converges should be independent of the initial condition. Therefore this result can 
be generalized to any initial conditions. Note that if our argument on the dominance of the 
effective potential is valid, this is an exact result in the infinite N T limit. 

This interpretation of the effective potential will enable us to determine the stationary 
distribution to which the complex Langevin equation converges. 

C. Effective Potential in One Loop 

In this section, we will calculate effective potential to first order in Q in the infinite 
Langevin time limit. As we showed in the previous section, in this limit the space-time and 
Langevin time average of the field must equal to the minimum of the effective potential. 

Calculation of the effective potential for a complex Langevin simulation that runs for a 
finite time is not easy. The calculation requires, as we will see, the evaluation of a 2N T N d - 
by-2N T N d matrix determinant and to the authors, there does not seem to be a way of 
simplifying this calculation. Instead we will do our calculation in the infinite Langevin time 
limit, but we will take this limit in a different way than what we have been considering until 
now. Until now, we let the complex Langevin simulation to start at k — with a fixed 
initial condition and run until a finite time N T . Then we took the limit N T — > oo. Instead, 
we now start our simulation at time k = —N T with a fixed initial condition and run it until 
N T . Then we let N T go to infinity. In the large N T limit, due to ergodicity and Langevin 
time independence of the complex Langevin equation, the Langevin time averages of these 
two processes will converge to the same stationary distribution. Another way to see the 
equivalence is by shifting the latter process by N T steps in Langevin time and thus reaching 
a process of the former type. Note that this changing of limiting process does not effect 
the Schwinger-Dyson equation related arguments. Since stationary distributions sampled 

4 We assume no degeneracy among the minima. 
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by the complex Langevin equation for both processes are the same, we can calculate the 
effective potential for the latter process, and its minimum will match the minimum of the 
former one. Matching of minima is sufficient for our purposes, but it is very likely true that 
the effective actions themselves match also. From now on, we work in the limit N T — > oo 
where the limit is taken in the latter way. 

We start by expanding the integrand of the generating functional for dynamical correla- 
tion functions around constant field configurations <p R and (p 1 . First, we choose the constant 
source fields as 



JRJ = 1 81$*, 0'] 



Using these values of the source fields, we get 



(35) 



expW[J H , J ] = exp 



a 



a — 1 



xc/"[#]exp -^E E 0^)M nKAm7B ^( 7 ) + O(0 3 ) 

L n,m k,7 A,B=R,I 



(36) 



Here the k sums go from — oo to oo. The matrix M is given by 



M, 



nnR,rri'yR 



dF 



R 



pa 



M, 



d(p R 

a \ d<f) R J ' af d<fi R2 ol - 1 \d(f) R J oi - 1 d(j) R2 

OF 1 



i >f' ; \ 2 i ~ R d 2 P R i //;/• ' V 2 i - ;!?/// 

' + —F — - — + f — ) + -F 



2e(a — 1) z — ' z(ct — 1) V 



■</ \ 1 ,i) 2 F' 1 f dF R Y 1 - r,d 2 F R 
+ -F 1 —^- + - — ^ + — .F 



M nKl 

,nnl 



1 IT) T) \ f)f?R 



a 



a — 1 



+ 



1 ^ 2 F' 

+ tF —r- — 4 



R^I 

1 0F fi 0F R 



1 * p 2 F R 

a 00«00' a-l' 00^0^ ' a(l-a)d(f) R 8(f) 1 



(37) 
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where 



k 

When deriving these equations, we made use of the Cauchy-Riemann equations. They are 
evaluated at the constant field configurations (f) R and (f) 1 . Now, using the definition of the 
effective potential, it is easy to see that 



F R (<I> R , 6 1 ) 



1 e 



F 1 (6 R M) 



In det M 



V(0 fi ,0 7 ) = -^ - + 7^^-n TT^+ n^ + 0(0). (39) 



n 4a n 4(a-l) 2£ n £ K 

Before continuing our calculation, we stop to discuss this result. As expected, the leading 
order term sets the minima of the effective potential at the classical fixed points, which are 
saddle points of the action S. This is just saying that if there were no noise (and the 
simulation did not diverge), the fields will evolve to a fixed point and stay there for ever. 
Averages over space-time and Langevin time in the infinite Langevin time limit will be set by 
the classical fixed points. At the classical level, all saddle points of S are degenerate, meaning 
they are all minima of the effective potential. Adding noise will remove the degeneracy 
through higher terms in the Q expansion. We will calculate the determinant term below 
and show in a specific example how it removes the degeneracy. 

On a different line of argument, we showed above that if the complex Langevin equation 
converges to a stationary distribution, this distribution will be a solution of the Schwinger- 
Dyson equation for the theory with complex action S. Moreover, we showed that Q counts 
the number of derivatives (loops) for this theory. Now we merge this result with the discus- 
sion of the previous paragraph as follows: 

Assuming ergodicity, averages over the stationary distribution is equivalent to averages 
over Langevin time in the infinite Langevin time limit. Since we set periodic boundary 
conditions on the space-time lattice, the resulting stationary distribution will have translation 
invariance in space-time. In particular, the Langevin time averages for the field variable 
will be constant in space-time. Therefore, we can say that Langevin time average of the 
field variable <p is equal to the Langevin time and space-time average of <p. Based on our 
results above, we know that the latter is equal to the minimum of the effective action for 
the complex Langevin equation which is given in equation (39) to th order in Q. Since Q 
simultaneously counts loop order for the stationary distribution, which is a solution to the 
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Schwinger-Dyson equations, we conclude that the complex Langevin equation converges to 
a stationary distribution that is quantized around a saddle point of the complex action S. 
This saddle point is given by the minimum of the effective action of the complex Langevin 
equation. 

The last statement needs more explanation. First of all, by quantization of a complex 
action S we mean not the Feynman path integral quantization but a Schwinger action princi- 
ple quantization. From Schwinger action principle, one derives Schwinger-Dyson equations, 
which have multiple solutions [10, 17]. 5 A solution of Schwinger-Dyson equations defines a 
generating functional for a quantization of S. This is consistent with equation (16), where we 
showed that in zero-dimensions solutions of Schwinger-Dyson equations can be represented 
by linear combinations of complexified path integrals. This result trivially generalizes to a 
lattice [4, 17]. 6 Complex Langevin simulations have only one stationary distribution, the 
question is to find which quantization of S this distribution gives. Our argument shows that 
the stationary distribution will contain contributions only from one of the saddle points of 
S. 7 One way of making this statement precise is through the results of [17]. There it was 
shown that for zero- dimensional theories, loop expansions of Schwinger-Dyson equations can 
be Borel resummed to give path integrals of type (13) with the path T being equivalent to a 
steepest descent path passing through a dominant saddle point of S. It is these generating 
functional for zero- dimensions, and their generalizations to higher dimensions that we argue 
to be the stationary distributions of the complex Langevin equations. Examples of this will 
be given in the next section. 

Now we go back to the determinant calculation. Calculation of the determinant of matrix 
M is most easily done in the Fourier basis. The discrete time Fourier transform and the 

5 Other differential equations also follow from Schwinger action principle, but it was shown in [17] that a 
study of Schwinger-Dyson equations is sufficient since in the continuum limit solutions of Schwinger-Dyson 
equations satisfy the remaining differential equations. 

6 To our knowledge generalizations of generating functionals with complex integration paths are not yet 
constructed on the continuum. Complex Langevin equations can be of help in defining these continuum 
generating functionals. 

7 Here we assumed that the th order term of the equation (39) breaks the degeneracy among the minima 
completely. If not, one may need to go to higher orders. If the degeneracy is not broken at all, it may 
be that the stationary distribution contains contributions from all degenerate minima. The generating 
functional for this distribution may be defined by a linear combination of the generating functionals that 
contain only contributions from one of the degenerate minima. See [17]. 
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inverse transform are given by 



(27r)( d+1 )/ 2 
1 



K 

P'K 

dw / dli . . 



-in-l 



dl d ' 



( 2vr )(rf+l)/2 J ^ 

The nonzero elements of the matrix M in this new basis are: 8 
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+ - 



+ 



dF 



Mwu/i =5(1 - I') 5(w - w' 



a «90«2 a -I \d(j) R 
1 
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dF 1 
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My 



n - r d(p 12 a \ dcf) 1 ) 'a' d(f) 12 

1 dF R ( V loj 



Mv«rm*i = 8(1 - S(w - w') 



V 



+ 



1 _ pi d 2 Fi 



2 dcj) 1 \ a a - 1 

1 dF R dF R 



\ P n d 2 F R 

a dj) R d& a-V d4> R d^ ' a(l - a) d(j> R <90 7 



Here, 



V lu} = e l 



k=i 



(40) 



(42) 



(43) 



It is now easy to calculate the logarithm of the determinant of this block- diagonal matrix 
where each block is a 2-by-2 matrix. 

det m 1 r dw r dh r dl d 



r dw r dh r di d , ^ 

/ — / 7; — • • • / 7— In [MiuRIurMiuI IuI — l-MiuRjul] J 
J-n 27T J_ n 27T J_ n ZTX 



2 En £ K 2 J-* 27r J-« ^ 



(44) 



recall the Fourier series representation of the S function: 



(41) 
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Using these results, the final expression for effective action in one loop is 



lim V(0 fi ,0 7 ) 

N T ^oo 



F R $ R , 7 ) 



4a 



+ 



F 1 (<p R , (p 1 ) 



n 



4(a - 1) 



+ 



i r dw r dh 



2ir 



2ir 



did 
2tt 



In (Mi^nM^^j - |M^ W | 2 ) + 0{Q). (45) 



In the next section, we will put equation (45) into use for specific examples. 



IV. AN EXAMPLE: ZERO-DIMENSIONAL CUBIC THEORY 



In this section, we will study zero-dimensional 3 theory as an example. The action is 
given by, 

S = -i4> - % -<p\ (46) 
while the Schwinger-Dyson equation is 

.d 2 Z 



— %■ 



df 



iZ(j)=jZ(j). 



(47) 



(48) 



We can construct the path integral solutions to this equation as: 

= j r d<f) exp {i<f) + |0 3 + j(j>) 
J T d(f) exp [i(f) + |0 3 } 

where T is a path that starts and ends in one of the wedges shown in Figure 1. Since 
the Schwinger-Dyson equation is second order, only two of these solutions are independent. 
One needs two boundary/initial conditions to fix a solution. One condition is given by the 
normalization Z(0) = 1. Setting another condition will fix the solution completely. 

Now we turn to the complex Langevin dynamics defined by the action (46). The dis- 
cretized complex Langevin equation is 



(f) R ( K + 1) = (f) R ( K ) - 2^t<P r (k)<P\k) + v/Att? 



Ri 



Noise correlators are 



At 



(49) 



(fi R (K)fi R ( K ')) = 2Qa5 KK/ 
(f, 1 ^) f l I {n')) = 2n{a-l)8 KK , 



(50) 



18 



FIG. 1: The region denned by sin(3#) > 0, where 9 is the argument of <j) is shaded. Each wedge is 
7r/3 radians wide. Any path starting and ending at the infinity within this area corresponds to a 
particular solution of ultralocal 3 theory. 



Classical fixed points of the complex Langevin equation is given by <p = i and <p = —i. As 
mentioned before, both fixed points are minima for the effective potential in the leading 
order. The next order for = % is given by 

1 



and for 



iln- 

2 4(Ar) z «(a-l) 
—i it is 

1 



r dw 



In 



4 (At) 2 + 4(1 - 2Ar) sin 2 - 



(51) 



2 "A_\2 



4 (At) a(a- 1) 



+ 



* dw , 
— In 



4(Ar)" + 4(l + 2Ar) sin 2 - 



(52) 



The integral in the above expression are convergent and <fi = i classical fixed point leads 
to a smaller term. Based on our discussion in the previous sections, we conclude that the 
complex langevin equation will converge a stationary distribution that is quantized around 
<p = i. To be more precise, we mean that the correlators of this theory are given by 



J r d(j) 4> k exp |i0 



(53) 



f T d(j) exp {i(f) + |0 3 } 

where T is the steepest descent path that passes from = i. This steepest descent path is 
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easy to identify is one changes variables to = —itp. Then 

_ H) fc / r , #^exp{^-^ 3 } 

J r ,#exp{^-^3} > ^ 

and now the saddle point that we want to quantize around is at ip = — 1. In this form, T' 
is recognized as the steepest descent path that defines an Airy function evaulated at 1, i.e. 
Ai(l). We refer the reader to [18] and [19] for details of this analysis. In particular, 

/ v Ai '(l) / X 

(0) = -< Ai(iy = 1Am > ( 55 ) 

which was calculated by MATLAB's built in Airy functions to 3 decimal places. In this 
equation, ' refers to a derivative with respect to the argument. Now, going back to equation 
(47), one concludes that 

<0 2 > = -l, <0 3 > = i- ((f)) = -0.176 ? , .... (56) 

((f)) fixes all the other correlators through algebraic relations. This is expected since setting 
((f)) to some number means imposing a second condition, aside from normalization, on the 
solutions of the Schwinger-Dyson equation (47). This fixes the solution completely. 

We expect a long time complex Langevin simulation to give us these correlators. We did 
100 complex Langevin simulations that ran from r = to r = 10000. The step size was 
At = 1 x 10~ 3 and a was chosen to be 1.001. Each simulation started from a random initial 
condition. The first three correlators calculated were 

(<fy) CL = 0.001(5) + 1.176(l)i, <0 2 > CL = -0.9999(5) + 0.003(13)i, 

((f) 3 ) CL = -0.006(21) - 0.176(7)i (57) 

Here the Langevin time averages are averaged over 100 simulations and the error bars stand 
for the standard deviation over 100 simulations. These numbers are in good agreement with 
theoretical results, as expected. 

We note that one could expect the (f> = % classical fixed point to dominate over (f> = —i 
classical fixed point by noticing that the former is a stable fixed point for the classical process, 
whereas the latter is not. Adding noise eventually makes the simulation sample around the 
stable fixed point, leading to a quantization around the stable fixed point. This parallels the 
results of [5, 6]. Figure 2 shows a set of sample points from the complex Langevin equation 
with classical flow on the background. The simulation is seen to sample around the stable 
fixed point (f> = i. 
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Re * 

FIG. 2: 1000 sample points (marked by +'s) from a complex Langevin simulation of action (46) 
with classical flow plotted in the background. 

V. CONCLUSION AND DISCUSSION 

In this paper, we attempted a new approach to the question of identifying the stationary 
distribution sampled by a complex Langevin equation. We worked in a lattice setting. The 
continuum limit and renormalizability issues are not trivial and deserves another study of 
its own. 

Our idea is based on a two step process. First, by solving a Schwinger-Dyson equation 
[4] one finds a possible set of stationary distributions. These are complexified path integrals 
[10, 17]. To identify the particular stationary distribution sampled by the complex Langevin 
equation, one needs to use information from the dynamics of the complex Langevin process. 
All the information about the dynamics of a complex Langevin process is in its generating 
functional. Our second step, whose details are given in this paper, is to use this generating 
functional and define an effective potential based on it. We argued that this effective po- 
tential governs the probability distribution of space-time and Langevin time averages of the 
simulated field variable. In particular, in the infinite Langevin time limit, the field average 
must equal to the minimum of the effective potential. We calculated this effective potential, 
expanding it in a parameter that multiplies the noise in the complex Langevin equation. In 
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the lowest order, the minima are given by the saddle points of the action S. The next order 
removes the degeneracy. The same parameter is shown to count loops in a loop expansion 
of the stationary distributions (that satisfy Schwinger- Dyson equations). Therefore, we con- 
cluded that the complex Langevin equation must be sampling a stationary distribution that 
describes a field theory that is quantized around the minimum of the complex Langevin 
equation effective action. 

As a final note, we want to mention that this construction can be used as a tool for 
studying symmetry breaking effects in stochastic quantization for complex actions [20]. 
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